Field evolution of the magnetic structures in Er2Ti207 through the critical point. 
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We have measured neutron diffraction patterns in a single crystal sample of the pyrochlore com- 
pound Er2Ti207 in the antiferromagnetic phase (T=0.3 K), as a function of the magnetic field, up to 
6T, applied along the [110] direction. We determine all the characteristics of the magnetic structure 
throughout the quantum critical point at Hc—2T. As a main result, all Er moments align along the 
field at He and their values reach a minimum. Using a four-sublattice self-consistent calculation, we 
show that the evolution of the magnetic structure and the value of the critical field are rather well 
reproduced using the same anisotropic exchange tensor as that accounting for the local paramag- 
netic susceptibility. In contrast, an isotropic exchange tensor does not match the moment variations 
through the critical point. The model also accounts semi- quantitatively for other experimental data 
previously measured, such as the field dependence of the heat capacity, energy of the dispersionless 
inelastic modes and transition temperature. 



I. INTRODUCTION 

Geometrical magnetic frustration allows one to obtain 
materials with tunable properties. It yields a large land- 
scape of possible magnetic ground states, due to the in- 
ability of the system to choose a unique spin configura- 
tion which minimizes the energy. Rare earth pyrochlores 
R2Ti2 07, where the R magnetic moments reside on the 
summits of corner sharing tetrahedra, are model systems 
to study such effects. Here, geometrical frustration does 
not arise from the competition of magnetic interactions, 
but rather emerges in the context of a highly symmet- 
rical structure, from the subtle interplay of three main 
energy terms: the single ion crystal field anisotropy, the 
exchange interaction and the magnetic dipolar coupling. 
In R2Ti207, the trigonal symmetry of the crystal field at 
the R site comes from the oxygen environment. At low 
temperature, this yields two generic behaviors, Ising-like 
(Ho, Dy, Tb) or XY-like (Er, Yb), depending on whether 
the <111> axes are easy or hard anisotropy axes for the 
magnetic moments. The final selection of a magnetic 
state within the ground state manifold is determined by 
the nature, lengthscale and sign of the magnetic inter- 
action, by perturbation energy terms, or by "order by 
disorder" processe^i^. 

The spin ice compounds Ho2Ti207 and Dy2Ti207, 
combining a strong Ising anisotropy with an effective 
ferromagnetic exchange, show exotic short range or- 
ders with macroscopic entropy^ and peculiar excitations, 
where the R dipolar moments fractionalize into mag- 
netic monopoles^. Here we focus on Er2Ti2 07, a com- 
pound with "reversed" behavior, namely with planar XY 
anisotropy and antiferromagnetic (AF) interactions, and 
showing magnetic order below 1.2 K^. 

Er2Ti2 07 was proposed to realize a model type XY 



antiferromagnet, for which theory predicts a fluctuation 
induced symmetry breaking, leading to magnetic long 
range ordering^. Indeed, for the crystal field ground 
Kramers doublet of the Er^+ ion, the threefold symme- 
try [111] axis is a hard magnetic axis, the easy plane 
being the local (111) plane. Below Tn=1.2K, the an- 
tiferromagnetic structure has a k=0 propagation vector 
and it is defined by the basis vectors ?/^2 , which transform 
according to the irreducible representation Fs, following 
the Kovalev notations used in Refs.6 and 7. It consists 
of six equally populated domains, as shown by spherical 
neutron polarimetry^. The selection of this particular 
state among the possible basis states of the k=0 mani- 
fold is still subject to discussion. Surprisingly, it is the 
only non-coplanar structure among all others, whereas 
order by disorder processes usually select coplanar or 
collinear orders^. It also differs from the the so-called 
Palmer- Chalker state^, predicted to be the ground state 
in the presence of AF isotropic exchange and dipolar in- 
teractions. An energetic selection of the ordered states 
was recently proposed^^, arising from sixth order terms 
in the crystal field Hamiltonian, together with possible 
anisotropic exchange, dipolar and Dzyaloshinskii inter- 
actions. An anisotropic molecular field tensor with an- 
tiferromagnetic components, reinforcing the crystal field 
planar anisotropy, has indeed been shown to be present in 
Er2Ti207 using polarized neutron diffraction, by measur- 
ing the local susceptibility in the paramagnetic phase^i. 

The magnetic structure near T=0 has been shown to 
coexist with spin fluctuations by /iSR spectroscopy^^: a 
non-vanishing spin dynamics washes out the precession 
signal usually observed by muons in ordered magnets. 
This is confirmed by the presence of soft collective exci- 
tations probed by inelastic neutron scattering^^. In this 
latter work, it is shown that application of a magnetic 
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field (along [110]) decreases the Neel temperature, result- 
ing in a zero temperature phase transition at a critical 
field He — 1.8 T. The magnetic ordered state is suggested 
to transform into some kind of spin liquid state above 
the critical field, through a second order quantum phase 
transition driven by spin fiuctuations^^^^^. The field evo- 
lution of the magnetic structure was qualitatively under- 
stood as a smooth deformation from the zero field ip2 
configuration to a field configuration above He where the 
moments are aligned with or close to the field direction. 
In these studies, however, the spin configurations could 
not be characterised in detail, considering the limited in- 
formation provided by the cold neutron measurements. 

In this work, we present detailed in-field neutron 
diffraction experiments in the AF phase of Er2Ti207 us- 
ing hot neutrons. This allows us to follow the evolution 
of the ground state induced by a magnetic field applied 
along [110], through the quantum critical point and up 
to a field of 6T, and to clarify both questions of the 
domain structure and of the microscopic local spin struc- 
ture within a tetrahedron. We show that the Er moment 
magnitude presents a minimum around He, which can 
be considered as the spin-fiip field for the [110] direc- 
tion. Up to 6T, the Er moments do not recover a fully 
collinear structure because [110] is not a principal direc- 
tion for the local ^-tensor of half the Er sites. We apply a 
self-consistent calculation, in the molecular field approx- 
imation, to try and explain this evolution quantitatively; 
our model takes into account the crystal field interaction 
previously studied in Refs.6 and 11, together with an 
anisotropic molecular field tensor, arising from exchange 
couplings with the 6 nearest neighbours of a given ion. 
We show that we can account for the field evolution of the 
magnetic structure in the AF phase (except in a limited 
field range below He) by the same molecular field ten- 
sor which accounts for the thermal variation of the local 
susceptibility in the paramagnetic phase^^. We also ap- 
ply our model to compute the lowest excitation energies 
at 0.05 K, the heat capacity in the paramagnetic phase, 
and the (i^, Tn) phase diagram; a reasonable agreement 
is obtained with the experimental data of RefstG, and.l3i. 



II. EXPERIMENTAL 

A single crystal of Er2Ti2 07 was grown by the fioating- 
zone technique, using a mirror furnace. It was placed at 
the bottom of a dilution inset, inside a superconduct- 
ing coil. The vertical axis of the magnetic field was 
aligned with the [110] axis, with a slight misorientation 
discussed below. Neutron diffraction measurements were 
performed at the ORPHEE reactor of the Laboratoire 
Leon Brillouin, on the Super-6T2 spectrometer in the 
unpolarized neutron version, with an incident neutron 
wavelength An=0.9 A. The nuclear structure was charac- 
terized by zero field neutron diffraction at two tempera- 
tures above Tn, 100 K and 5K, allowing the lattice con- 
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FIG. 1. Integrated intensities versus the field H applied 
along [110] of some typical Bragg peaks measured at 0.3 K in 
ErsTisOr. 



stant, positional parameters, occupancy factors, isotropic 
temperature factors and extinction parameters to be re- 
fined, within the space group Fd3m. To determine the 
magnetic structures, about 300 Bragg peaks were col- 
lected at 0.3 K, in zero field and for 10 field values in the 
range - 6 T. 



III. EVOLUTION OF THE MAGNETIC 
STRUCTURE 

Figure [T] shows the field dependence of the integrated 
intensities of some particular Bragg peaks. The (/i, — /i, /) 
peaks, situated in the horizontal plane, have a field de- 
pendence similar to that previously measured^^. The lift- 
ing counter geometry of the diffractometer allowed us to 
measure the equivalent out-of-plane (h^h^l) refiections 
plotted in the figure, which show different field depen- 
dencies. These allow several regions of interest to be 
distinguished. The low field region H < 0.1 T shows 
opposite variations of the (220) and (022)/(202) Bragg 
peaks, which are clearly assigned to the reorientation of 
the magnetic domains by the applied field. At higher 
fields, a single domain is stabilized, and the variations of 
the Bragg peaks arise from moment reorientations inside 
a tetrahedron. The (002) and (111) Bragg peaks show an 
extremum at 1.5 T, somehow below the critical field. The 
evolution of the magnetic peaks with the field is well ac- 
counted for by the magnetic refinements described below. 

Magnetic refinements were performed with the pro- 
gram Fullprof^^. In zero field, the magnetic structure is a 
k=0 structure, which means that the 4 tetrahedra in the 
cubic unit cell have the same moment orientations. These 
orientations are those corresponding to the state, with 
moments in the easy planes along <211> axes. The mag- 
netic moment at 0.3 K was refined as m=3. 25(20) jub per 
Er ion. The domain populations were refined, yielding 
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FIG. 2. [Color on line] Zero- field configuration in the 1^2 state 
in the (110) (left) and (110) (right) plane. The [110] incipient 
field direction is indicated by a green arrow. The local ternary 
axis (hard magnetic axis) for the ai (resp. 0^2) moments 
is [111] (resp. [Ill]); for the /^a (resp. (S^) moments, it is 
[111] (resp. [111]). This configuration corresponds to that 
of the single domain occurring above 0.1 T. The upper (resp. 
lower) hemisphere referred to in the text is that containing 
the positive (resp. negative) part of the [001] axis. 



three equally populated magnetic domains (together with 
the opposite domains giving the same contribution), as 
also found by polarized neutron measurements^. 

Above 0.1 T, good refinements (with typical agreement 
factors Ri?=5%) were obtained by considering a single 
magnetic domain. The evolution of the magnetic struc- 
ture with the field was determined by refining the in- 
tegrated magnetic intensities, with moment values and 
angles as parameters. For the [110] field direction, the 
Er sites split into two sets, or chains: the a sites, with a 
ternary axis at an angle 6 = arcsin(l/>/3) 35.3° from 
H, and the /3 sites with their ternary axis perpendicu- 
lar to the field. In the ?/^2 state, the a sites split further 
into ai and a2, with moments which are not symmetrical 
with respect to the field direction. We also split here the 
/3 moments into Ps and (see Figl2]). In our experimen- 
tal setup, there is a slight misorientation of the applied 
field with respect to the [110] direction: the polar angle 
of the field is 93.8° (instead of 90°) and its azimuthal 
angle cp is 45.8° (instead of 45°). In the calculation de- 
scribed in the following section, we will neglect the small 
azimuthal misorientation and assume (f = 45°, but we 
take into account the exact value. In the refinements, 
in order to reduce the number of fitting parameters, we 
assumed equal magnitudes for the two /3 moments, which 
was a posteriori justified by the calculation. The pre- 
ferred domain orientation (above a field of about 0.1 T) 
is such that the a moments lie in the upper hemisphere. 

The moment orientations in a given tetrahedron are 
sketched for selected fields in Figl3l All the moments are 
seen to rotate towards the field direction (see c or d in 
FigH]) up to a field of about 2 T, which can be identified 
with the quantum critical field He of Ref. 13. For this field 
value, all the moments are practically aligned along the 
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FIG. 3. [Color online] Evolution of the spin structure inside 
a tetrahedron for selected fields. The flat disks visualize the 
easy (111) planes. The a moments are shown by red arrows 
and the (3 moments by blue arrows, with lengths proportional 
to the moment values. 



field, so that He can be considered as the spin-flip field 
of the Er2Ti207 magnetic structure for the [110] field 
direction. On further increase of H beyond 2 T, the /3 
moments remain along the field whereas the a moments 
tend towards an asymptotic orientation, symmetric with 
respect to H. 

The data points in Fig H] represent the evolution of the 
moment values and orientations at the 4 Er sites as the 
field increases. Pronounced anomalies are seen in the 
quantum critical region. As a major effect, all moment 
values show a minimum at Hc=2T (FigHJa). This mini- 
mum is much more pronounced for the a moments, which 
get the closest to their hard local ternary axis for this field 
value. The ai moment decreases down to 1 /ie whereas 
the a2 moment reaches 2 /iB- As to the P moments, which 
remain in their easy plane while rotating, they show a 
much shallower minimum of about 3 /iB • 

The evolution of the angles of the moments with the 
field is shown in FigsH] c and d. At 2 T, all the angles 
are close to zero. The "flip" of the ai moment is re- 
flected by the change of sign of its angle with H, when 
it enters the lower hemisphere (see Figl2j) above 2 T. By 
contrast, the a2 moment approaches the field direction 
at 2T, then tilts away, always remaining in the upper 
hemisphere. Above 2T, the two a moments tend to- 
wards an orientation at 20-25° on either side of the field. 
The P moments show a much smoother field variation: 
they progressively reorient towards the field and remain 
aligned along H from 2T upwards. The angles of the 
moments with their local ternary axis (referred to as OZ 
in FigH] b) confirm that the P moments always remain 
in their easy plane, while the a moments approach their 
hard axis in the critical region. 

Our model calculation, to be described in the next sec- 
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FIG. 4. [Color on line] Variation with the field, applied along [110], of the moment magnitude (a) and orientation with respect 
to the local ternary axis (b) and to the field direction (c and d), for the 4 Er sites in Er2Ti207 at 0.3 K. Our convention for 
the angle with the field direction is that it is negative if the moment lies in the lower hemisphere. Solid lines in (a)-(c) are self- 
consistent calculations with = —0.51 Tf/HB and Ay = — 0.06T//is; Solid lines in (d) are calculations with Ax = — 0.54T//is, 
Xy = — 0.48T//XS and X\\ = — 0.06T//is; dashed lines are interpolated ones. 



tion, accounts for most features of this evolution. 



IV. MODEL CALCULATION 

We implemented a model intended to reproduce the 
experimental data concerning both the thermal variation 
of the local magnetic susceptibility in the paramagnetic 
phase, measured in ReflTTI, and the evolution with field of 
the magnetic structure at 0.3 K, measured in the present 
work. This model performs mean field self-consistent cal- 
culations and uses as ingredients the cry stal field param- 
eters of Er^+ in Er2Ti207 as in Refflll. and anisotropic 
two- ion exchange of the type: 

y^ex = —J\\SizS2Z — J^{SixS2X + SiyS2y)^ (1) 

where J\\ and are the components of the exchange 
tensor in the local frame with axial symmetry. We shall 
also consider the effect of a slight exchange anisotropy 
within the easy plane. We consider each ion in a tetrahe- 
dron to be exchange-coupled to its 6 nearest neighbours, 
resulting in an anisotropic molecular field tensor A such 
that the molecular field acting on this ion writes: 



1~ ^ 
Hmoi = TrAy^nife, 



6 



(2) 



where the sum runs over the 6 nearest neighbors. The 
relationship between the components of the A and J ten- 



sors is: 



Ai = ^Ji ( 



1 



1 



(3) 



k=l 



A self-consistent treatment involving the 4 Er moments, 
each with its 3 components, is performed, the only pa- 
rameters being thus the 2 components of the A-tensor. 
The dipolar coupling is not included in the calculation. 

We apply first this model to the paramagnetic suscep- 
tibility and to the zero-field AF phase. In the param- 
agnetic phase, the X||(^) ^^^^ X±(^) data measured in 
Refilli are very well reproduced using the above model 
with the axially symmetric AF molecular field tensor: 
A|| = -0.06(3) T//iB and Ax = -0.51(4) T//iB (see 
FiglSj). These values are close to those obtained in RefUll 
using a simpler single sublattice model (Ay = —0.15, 

Ax = —0.45 T//iB)- The planar anisotropy of the A 
tensor is strong (Ax/Ay :^10) and reinforces that of the 
crystal field. In the AF phase in zero-field, the 4 Er 
moments in the ?/^2 ground state comply with the rule: 
mj=0, whence, for instance, m2 + ma + m4 = —mi. 
Therefore, the molecular field acting on ion i is, ac- 
cording to Equ.Q: HJ^^^ = — ^A^m^, since the spon- 
taneous moments lie in the easy plane. According to 
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FIG. 5. [Color on line] Low temperature region of the thermal 
variation of the longitudinal and transverse susceptibilities 
in Er2Ti207. The data are taken from ReflTTI and the lines 
are calculated using the self-consistent model described in the 
text. 



mean field theory, the transition temperature is given 
by: /cbTn = ^o^moh hence in the ip2 state of Er2Ti207: 
/cbTn = 1 1 Ax I 5 where the saturated moment mo and 
X_\_ must be obtained self-consistently in zero applied field 
and near OK. In order to obtain a transition tempera- 
ture of 1.2 K, one needs the value = — 0.435 T//iB, 
which is remarkably close to that derived in the para- 
magnetic phase. This coherence gives confidence in the 
applicability of the molecular field approximation in the 
paramagnetic and AF phases of Er2Ti207. The associ- 
ated transverse exchange integral is J7± = — 1.75K, and 
the T = molecular field H^oi = 0.5 T. The zero-field 
spontaneous moment is mo = 3.52 /ie, close to the ex- 
perimental value 3.25(20) /xb, but somewhat higher than 
3.00(05) /iB measured in Ref.6. 

We then used the self-consistent exchange and crys- 
tal field model described above to calculate the evolu- 
tion of the magnetic structure with increasing field. In 
a first step, we chose the two components of the axi- 
ally symmetric A tensor as determined from the fit of 
the paramagnetic susceptibility {X± = — 0.51T//iB and 
A|| = — 0.06T//iB)- The results are displayed in FigsHJa, 
b and c as solid lines. The calculation does not exactly 
match the data, but it captures the main trends of the 
rotation of the Er moments as the field is increased. Re- 
markably again, using the A tensor derived in the param- 
agnetic phase, the value of the critical or "spin-flip" field 
He ^2T is correctly reproduced. This value strongly 
depends on A^ and very little on Ay. For a simple AF 
structure with moments perpendicular to the field, the 
spin-flip field in the presence of anisotropic g and A ten- 
sors writes: 



i^sf = ^^||/iB |A|| + (^)^Ax|. 

^ 3\\ 



(4) 




Since its moment configuration is more complicated than 
a two-sublattice AF structure, this expression cannot 
be directly applied to Er2Ti207, but it accounts for 



FIG. 6. [Color on line] Evolution of the Er magnetic moments 
with the magnetic field at 0.3K, as calculated from the self- 
consistent calculation with with Xx = — 0.54T//is, Xy — 
— 0.48T//iB and X\\ = — 0.06T//iB- Arrows are proportional 
to the moment values. Each field step corresponds to 0.48 T. 



the strong dependence of i^sf on A^ since g^/g\\ >1 in 
Er2Ti207. 

However, as seen in FigHJc, there is a discrepancy con- 
cerning the rotation of the j3 moments: the calculation 
yields different variations of the angle with the field for 
the two sites, whereas the data present identical varia- 
tions within experimental uncertainties. This is due to 
a moment reorientation occurring in the calculation at a 
very small field (about 0.03 T): the moments jump from 
the ip2 state, with moments along the local OX axis of 
the <211> family, to a configuration where each moment 
lies at about 38° from OX, which probably corresponds 
to the (small) potential well created by the crystal field 
interaction. Since this reorientation is not observed in 
the data, one has to devise a mechanism which would re- 
inforce the potential well along the local OX axis, at least 
at low field. For this purpose, in a second step, we lifted 
the in-plane degeneracy of the molecular field tensor by 
about 5%, through the introduction a slightly higher (ab- 
solute) value for Ax than for Ay. The result of the cal- 
culation, with Ax = -0.54T//iB and Ay = -0.48T//iB, 
is shown for the angle ^H=(m,H) in FigH] d (the other 
calculated quantities remaining practically unchanged). 
The agreement is better for the P moments, i.e. the ip2 
configuration evolves smoothly with increasing field, up 
to a field of 1.2 T. For 1.2 < < 2T, the calculation 
presents an anomalous moment jump, still unexplained. 
Therefore, in FigH] d, we do not consider the calculated 
curves in this field interval and we have drawn dashed 
lines interpolated between the lower and higher field re- 
gions where the model reproduces the data reasonably 
well. In Figini we show the evolutions of the four mo- 
ments versus the magnetic field, calculated within this 
last model. 

Finally, the high field configuration is the following: 
the 13 moments lie along the field and the ai and 0^2 mo- 
ments are on either side of the field, at an angle of about 
20°. This can be understood since the [110] direction lies 
within the easy plane for the p moments, whereas it is 
not a principal axis for the a moments^^. 

It is worth emphasizing here that the data cannot be 
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FIG. 7. [Color on line] Calculated lowest excitation single 
ion energies for the 4 sites in Er2Ti2 07, in the AF ordered 
phase at 0.05 K as a function of the applied magnetic field. 
The square dots are the (estimated) experimental values from 
Ref 13, for fields where a clearly non-dispersive level is ob- 
served in the inelastic neutron spectra. Insert: calculated 
thermal variation of the heat capacity for i7=1.5, 3 and 7T, 
obtained from the single ion energies. 

reproduced with an isotropic molecular field constant A = 
—0.51 T//iB- In this case, the calculation yields moments 
which remain in their easy plane up to the critical field 
of 2 T, and this field is not a spin-flip field since the a 
moments do not align along its direction. The a moment 
magnitude then shows no minimum at 2 T. Therefore, 
our calculation, within its limits of validity, demonstrates 
the existence of an anisotropic molecular field tensor in 
ErsTisOr. 



V. DISCUSSION 

In this section, we show that the model described above 
also accounts relatively well for other experimental re- 
sults related to the spin excitations in Er2Ti207, as pre- 
viously measured inelastic neutron scattering and heat 
capacity. In Refs.6 and Ts", a dispersionless mode was 
observed at about 0.4 meV for i^=0, which further splits 
into higher energy modes when the field increases above 
the critical field. In first approximation, mostly valid in 
the middle of the Brillouin zone, such flat modes may be 
attributed to the splitting of the ground Kramers doublet 
under the coupled influence of the external and exchange 
fields. In this picture, the splitting into several modes 
induced by the field can be understood by the different 
molecular fields experienced by the a and /3 moments. 
In Figlll the excitation energies obtained by our model 
are compared with the energies of the flat modes deter- 
mined in Ref.13 (square dots in Fig|7j). Below 2T, the 
calculated energies are almost field independent around 
0.25 meV, and above this threshold field, they increase 
linearly with slightly different slopes. The data match 
reasonably well the calculated values, except at zero and 
low field, where the energy of the non-dispersive mode 




Temperature (K) 

FIG. 8. [Color on line] Calculated thermal variation of the 
ai Er moment in the presence of a magnetic field applied 
along [110], for different field values below 1.3 T, upper limit of 
validity of our model in the AF phase. Insert: (if, Tn) phase 
diagram obtained from the break-points of the ai-site m(T) 
curves. The dashed line represents the law H = Hc[l — {^)'^], 

N 

with ifc=1.85T and T^=1A8K. 



(0.4 meV) is significantly above the calculation. One ex- 
pects that our single ion mean field energy calculation 
holds at high field, where exchange is small with respect 
to the Zeeman energy, but fails at low field. Indeed, 
at low field, one should consider, for instance, the sys- 
tem of 4 exchange coupled Er moments on a tetrahedron 
and calculate the whole spin-wave spectrum, taking al- 
together crystal field, exchange and Zeeman interactions 
into account. 

In Ref.13., the variation with field of the heat ca- 
pacity in Er2Ti2 07 was attributed to a crossover to- 
wards a quantum high field paramagnetic state remi- 
niscent of that observed in the quantum critical mag- 
net LiHoF4. We calculated the in- field heat capacity 
Cp{T) in Er2Ti207 below 10 K, using the crystal field 
and Zeeman energies of all the levels (see insert in Fig|7]). 
Whereas this simple Schottky-type calculation cannot 
capture the anomaly at the phase transition for low fields, 
it reproduces correctly the main features of the high field 
heat capacity data of Ref.13. Especially, the smearing 
of the temperature dependence of Cp{T) with increasing 
field can be correctly reproduced in this simple approach, 
which does not involves quantum fluctuations explicitly, 
besides the crystal field interaction. 

Finally, we also applied our model to the calculation 
of the ordering temperature as a function of the applied 
field, i.e. to the (if, Tn) phase diagram. With the A- 
tensor determined above {Xx = —0.54, Xy = —0.48 and 
A|| = — 0.06 T//iB), we calculated the thermal variation of 
the Er moments for different field values in the tempera- 
ture range 0.2 - 2 K, in the field range below 1.3 T, where 
the model yields good agreement with experiment. As 
can be seen in FiglHl which displays the thermal variation 
of the ai moment magnitude, the ordering temperature, 
defined as the temperature where the moment variation 
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shows a jump, decreases with increasing field^^. This re- 
flects the competition between the Zeeman and exchange 
energies, which leads to the quantum critical point when 
the Zeeman coupling overwhelms exchange and prevents 
a spontaneous moment configuration to set in. The low 
field part of the (i^, Tn) phase diagram is shown in the 
insert of FiglHl Our calculation is in qualitative agree- 
ment with the diagram obtained from the heat capacity 
data of RefB The curve H = f{T^) is seen to follow 
a power law (see dashed line in the insert of Fig IS]), with 
an exponent n = 2. Extrapolation of this law to zero 
yields a critical field He = f{T^ = 0)=1.85T, close to 
the experimental value. 

VI. CONCLUSION 

Using neutron diffraction in the AF phase of the planar 
pyrochlore Er2Ti207 with a magnetic field applied along 
[110], we performed a quantitative study of the field evo- 
lution of the magnetic structure throughout the quantum 
critical point at He —2 T. For this critical field value, the 
AF magnetic structure has "fiipped" , i.e. all the Er mo- 
ments are aligned along the field and their magnitude 
reaches a minimum. The strong decrease of the moment 
magnitude at He refiects the spin wave damping and 
enhanced quasi-elastic scattering previously observed in 



neutron experiments. A four sublattice self-consistent 
calculation, taking into account exchange, Zeeman and 
crystal electric field interactions, accounts for most of 
the characteristics of the field induced magnetic struc- 
ture, from zero field up to well above He. It also explains 
semi-quantitatively the field dependence of the heat ca- 
pacity, of the dispersionless inelastic mode and of the 
transition temperature. The comparison between model 
and experiment brings out the strongly anisotropic ex- 
change interaction as a necessary ingredient to explain 
all the features of the Quantum Critical Point. The crit- 
ical field can then be defined as the field value for which 
the Zeeman energy overcomes the exchange energy, tak- 
ing the anisotropic exchange tensor into account. The 
anisotropic exchange, already outlined by local suscepti- 
bility measurements in the paramagnetic phase, should 
also infiuence the spin-wave excitation spectrum in the 
AF ordered phase. 
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